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Abstract 



> 

Competing magnetically ordered structures in polymerized orthorhombic 
^ ■ ACqo are studied. A mean-field theory for the equilibrium phases is de- 

t-H ! 

veloped using an Ising model and a classical Heisenberg model to describe 

On , 

the competition between inter- and intra-chain magnetic order in the solid. 
In the Ising model, the limiting commensurate one-dimensional and three- 
dimensional phases are separated by a commensurate three-sublattice state 

> 

^ ■ ami by I wo sectors containing higher-order comruonsurat.o phases, hor (hie 

■ Heisenberg model the quasi-lD phase is never the equilibrium state; instead 

the 3D commensurate phases exhibits a transition to a continuum of coplanar 
spiral magnetic phases. 
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Recent experiments have revealed that ACqq (where A=K, Rb) undergoes a first-order 
structural phase transition from a high-temperature face centered cubic (fee) structure to a 
low-temperature orthorhombic phase [|1[]. The structure of this new low-temperature phase 
is quite interesting since the lattice parameter along the polymer direction (c axis) is only 
9.1 A. This distance implies a separation between nearest-neighbor molecules along this axis 
as small as ~2 A, and has led to the suggestion that the fullerene molecules in this phase 
polymerize to form an ordered crystalline array of quasi-one-dimensional (quasi-lD) chains, 
resembling a lattice of aligned "necklaces" of fullerene molecules. Moreover, below 50 K this 
orthorhombic phase undergoes a second electronic phase transition, which Chauvet et al. 
argue arises from a spin- density- wave instability within the quasi-lD chains [0. 

More detailed electronic structure studies do not support this interpretation 0. The 
conduction band calculated for the orthorhombic phase shows substantial dispersion for the 
doped carriers both parallel and perpendicular to the polymer direction. Furthermore, by 
treating the effect of a short-range repulsive potential on this paramagnetic state within 
the random-phase approximation (RPA), one finds that the metallic phase of this structure 
favors magnetic fluctuations to a three-dimensionally ordered two-sublattice state, with spin 
polarizations alternating on the corner and body-centered sites of the pseudo-body-centered 
orthorhombic cell 0. Indeed, this spin structure is stable for quite modest values of the 
repulsion strength, and is likely associated with the 50 K transition observed experimentally. 

However, it is plausible that a dilation of the interchain galleries, such as might be 
achieved by alloying the alkali donors with larger spacer molecules, could lead to the one 
dimensional situation envisioned by Chauvet et al. If so, the transition of the magnetic 
structure from the three-dimensional ordered phase to the quasi-lD phase is of some impor- 
tance, and is the subject of this paper. We study the transition from three-dimensional to 
one-dimensional magnetic structures using both an Ising model and a classical Heisenberg 
model to describe the spin dynamics. Interestingly, as we show below, the passage from the 
three-dimensional ordered structure to the quasi-lD structure is not an easy one. In fact, 
at T = the quasi-lD modulated magnetic phase appears as an equilibrium phase in the 
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Ising model but not the Heisenberg model for this system. In the former case, the ordered 
quasi- ID magnetic structure is separated from the three-dimensional structure by an assort- 
ment of intervening modulated magnetic phases, and in the latter case it is preempted by a 
continuum of spiral magnetic phases. 

The structure of the pseudo-body-centered orthorhombic (pseudo-bco) phase is sketched 
in Fig. |](a). The measured lattice parameters (for RbCeo) are a=14.23, 6=10.11, c=9.14 
A ||. There are two candidate orientational structures for the fullerene molecules on this 
Bravais lattice; the measured x-ray diffraction intensities at present are unable to distinguish 
between a model with fixed orientations (space group Pmnn) and a model with binary 
orientational disorder (space group Immm) ||. For simplicity, we assume here a true bco 
Bravais lattice with a single orientation ||. The Brillouin zone for this bco cell is sketched 
in Fig. [l|(b). In our model, each Bravais lattice site of this structure is presumed to localize 
a single electron, with the coupling between the sites described by the Hamiltonian 



Heisenberg model (A = 1), with the spins treated classically. 

The sum in Eq. ([!]) runs over the first ten sites j surrounding any lattice site i. This 
sum thus includes all the interactions J\ between a reference site and its eight neighbors 
along the body diagonal of the bco cell, and the interaction J 2 to the two nearest neighbors 
along the polymer direction. Note that this spin Hamiltonian can be frustrated, since for 
the choices J\ > and J 2 > antiferromagnetic correlations along the body diagonal must 
compete with antiferromagnetic correlations along the polymer chain. For the actual RbCeo 
structure, our calculations within the RPA show that the system is in fact not frustrated 
0, and that both of the effective interactions favor the two-sublattice structure discussed in 
Ref. ||, that is, we have a situation with J\ > and J 2 < 0. More distant interactions within 
each of the chains also favor ferromagnetic intrachain order. However, one can imagine that 
that by reducing the size of the interchain hopping amplitudes one ultimately arrives in the 




(1) 



where Sf 



Sf ± iSf. We consider the limiting cases of the Ising model (A = 0) and the 
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regime with J2 >> J\ > 0, which would then favor the quasi-lD spin configuration originally 
proposed by Chauvet et al 0. In the following we label these two limiting commensurate 
reference states CVd (for the three-dimensional state) and C\d (for the one-dimensional 
state). Schematic representations of these spin configurations are shown in Fig. 0. 

The general structure of the phase diagram for these colinear states, shown in Fig. |3](a), 
can be understood by comparing the free energies of the competing colinear states, Cm 
and C3D, as we vary the interaction ratio J2/ J\- Within the Ising model, mean field theory 
yields a multicritical line at J2/J1 = 2, which separates the C^d phase from the Cm phase. 
Each of these phases is in turn separated from the high-temperature disordered phase by 
a second-order phase boundary, denoted by solid lines. There are an infinite number of 
coexisting phases along vertical dashed line, since the line locates the critical interaction 
ratio where the creation energy vanishes for a domain wall which slips the stacking sequence 
of the laterally ferromagnetically ordered layers by one lattice spacing along the polymer 
axis. 

The situation is more interesting if one allows for the possibility for a longer range 
modulation of the spin texture. Consider first the Ising limit, given by A = 0. The mean- 
field inverse susceptibility in the paramagnetic phase of the model has the form x _1 (?) T) = 
S%J(q) + kT, where S% = 1/4 and J(q) = J2j Jij exp(zg • Rij). For a given value of J%j Ji, we 
then seek the critical wave vector q c which minimizes J(q)- Examination of the evolution 
of q c reveals that this system exhibits four limiting instabilities from the paramagnetic 
phase: (a) For J2/J1 < 1, the dominant fluctuations are at the X point of the bco zone, 
(27r/a,0,0); below the critical temperature, this leads to the C^d spin structure, (b) For 
1 < J2/J1 < 2, the ordering wave vector evolves continuously from the X point towards 
the G point, (2n/a, 0, 27r/3c). (c) For J2/J1 = 2, the dominant fluctuation is obtained 
exactly at the G point. This leads to a commensurate three-sublattice antiferromagnetically 
ordered phase, C*3 S , which will be briefly discussed below, (d) For 2 < J2/J1, the modulated 
structure is again incommensurate, with the ordering wave vector evolving between the G 
point and the point, (2ir/a, 0, vr/c). The one-dimensional limit of this problem corresponds 
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to a condensation at any wavevector on the plane q z = ±7r/c. 

For < T < T c we can investigate the properties of the low-temperature ordered phases 
by solving self-consistently the linear stability relations o~j = <7 tanh(/3/ij), where the field 
is hi = — J2j Jiji&j)- The resulting phase diagram is shown in Fig. §|(b). As one expects in 
the frustrated Ising model, the degeneracy along the original multicritical line is broken by 
the entropic contributions to the free energy [Q, and at finite temperature one finds that a 
T = multicritical point branches into two sectors containing modulated incommensurate 
(actually, higher order commensurate) phases, I and and a locked commensurate three- 
sublattice phase, C^ s . This intermediate commensurate phase is a modulated phase with the 
amplitude o z {z) = a sm(2iiz/3c). It thus consists of lateral planes of fullerene molecules 
ferromagnetically aligned within a layer, with antiferromagnetic ordering across isolated 
pairs of neighboring layers, leaving the third layer in the cell unpolarized; this is shown 
schematically in Fig. 0. The regions / and I' are themselves further partitioned into a 
denser set of higher order commensurate phases. 

We now remove the scalar constraint on the spin degree of freedom, and treat the Heisen- 
berg vector-spin model of Eq. (|T]) with A = 1. The orientation of the spin on the i-th lattice 
site can be specified by an angle Qi locating a point on the surface of the unit sphere. Within 
the mean field theory, expanding around the disordered phase, we consider a normalized trial 
density matrix of the form p(Q) = (l/47r) (1 + J2 a PaKa(ty) with a = x,y,z labeling the 
three real normalized vector harmonics (K x , K y ,K z ) = y3/4ft (sin 6 cos 0, sin 6 sin 0, cos 6) . 
The thermal average of the a-th component is thus (S a ) = ^Jl/12n S Q p a . The mean-field 
free energy per site in this model has the form 

/ = ^TT7 Jij,a(Sia) (Sj a ) 
Zly ij, Q 

kT r 

+^rE/ dn lP (n t )\o gP (n t ). (2) 

i,a 

To study the equilibrium magnetic states of the system, we minimize this trial free energy 
with respect to the coefficients p a . For T << T c the expansion in K a is inadequate, and we 
replace it by an expansion in a set of gaussians uniformly distributed over the surface of the 
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unit sphere. 

To quadratic order in the p a , the free energy in this vector model takes the form: 

/ = ^5><«)*<^>P^ + k'). ( 3 ) 

The susceptibility thus diverges on the same locus as found above for the Ising model, and 
leads to the same high-temperature phase boundary as shown in Fig. |](b). The vector char- 
acter of the spin qualitatively changes the low-temperature behavior of the model, however. 
At T = the commensurate C^d phase remains stable only up to the critical point J2/J1 = 
1. This can be seen most directly by studying the long wavelength spin fluctuations around 
the ordered C^d state. To do this we adopt a semiclassical model, in which the equation of 
motion for a spin precessing on site % is: 

^ = S< x Kff, (4) 

where h^^ ia = —(dH/dS ia ), and the brackets denote a thermodynamic average. Eq. (f|) 
is then solved by linearizing the spin, S ~ (S) + er(£). The dependence of the equation 
of motion on the geometric structure of the parent spin state is expressed explicitly in 
Eq. (^). Expanding around the ordered C^d phase one finds that the competition between 
Ji and J2 renormalizes the spin-wave velocity. For q —>■ the spin waves have the dispersion 
lu = 8( Ji( Ji — J 2 ))^g 2 c , so that at J 2 / J\ = 1 the system is only marginally stable with respect 
to transverse fluctuations of the magnetization that are modulated along the c direction. 

For larger values of J2/J1, the system spontaneously develops a static transverse un- 
dulation of the ordered moment. This yields a helical phase modulated at a wave vector 
q that can be tuned continuously by varying the interaction ratio. Indeed, for J%j J\ > 1 
the equilibrium phase of the model is a coplanar spiral texture. This is seen by studying 
the zero-temperature internal energy per site, u. We assume that each spin is tipped with 
respect to the c axis by a polar angle 8, and that the transverse component of the magne- 
tization advances by a uniform pitch between neighboring planes normal to the polymer 
direction. The energy per site then is then 
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u = J2 cos 2 9 + J 2 sin 2 9 cos 20 + 4 Ji sin 2 # cos (p. (5) 

We see that this interaction energy is always minimized by a texture with 9 = it/2 and 
(ft = cos _1 (— J1/J2). The dependence of on J2/<A is shown in Fig. |j. One sees that for 
J2/J1 1 + the spiral phase smoothly connects to the ordered C%d phase (this corresponds 
to the choice 9 = n/2 and <fi = n). A second critical point occurs at J2/J1 = 2. Here <fi = 
27r/3 which describes the three-sublattice spiral phase. In fact, the geometric structure of 
the coplanar three-sublattice spin texture can already be deduced by a careful inspection 
of the sixth-order invariants in the Landau expansion of the free energy for a vector order 
parameter near the critical point connecting the three-sublattice phase to the paramagnetic 
phase. 

Remarkably, the energy of Eq. @ requires that for any finite J%j J\ (that is, for any 
nonzero J\) the ground state will always be a helical structure, albeit with a pitch which 
diverges as J2/J1 — » 00. This is confirmed by studying the spin- wave expansion for small 
fluctuations around a putatively ordered phase Cm at T = in the regime J%j J\ > 2. Here, 
global rotational invariance requires that the q = mode occurs at precisely u — 0, but 
for nonzero q near the zone center, we find that for all coupling strengths the transverse 
fluctuations around this phase are unstable. The width of this instability in momentum space 
collapses as J2/J1 diverges. Thus the commensurate structure C\r> is intrinsically unstable 
at T = for all finite coupling ratios. In Fig. f| one sees that the pitch asymptotically 
approaches the value tt/2 , which is its value for the Cm phase, only in the limit of large 
J2/J1. 

The phase diagram for the classical vector-spin model in the bco structure is shown in 
Fig. 0(c). The commensurate C 3 e> phase is stable for J2/J1 < 1, and the hatched regions to 
the right describe the continuum of coplanar spiral phases by which the C%d phase winds 
out, asymptotically approaching the Cm phase. 

The energy density of Eq. (^) requires that the helical phases are all coplanar spiral 
textures; thus nowhere in the phase diagram do we obtain a phase in which the spin texture 
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orders in three dimensions. However, this can be the case if we treat the effects of a crystal- 
field anisotropy. Orthorhombic symmetry allows an on-site potential of the form aS* + bSy + 
cS^. Competition between this on-site potential and the planar spiral textures favored by 
the quadratic interactions can then "tip" the spins out of the plane, yielding a stable three 
dimensionally ordered conical spin configuration. Of course, in the presence of an external 
magnetic field, a weak external field acts as a director which orients the plane of the spiral 
modulation normal to the applied field, with the spins canted uniformly along the external 
field direction. 

The analysis given above has been carried out for a model in which the spins are presumed 
to be localized on each lattice site of the structure. However, it is likely that the doped 
carriers of the polymerized solid are better described by the itinerant limit of this theory, 
along the lines of the theory discussed in Ref. 0. In fact, we note the experimental transition 
temperature to a magnetically ordered phase occurs at 50 K, while the Fermi energy for the 
doped carries is of order 3000 K, so that quantum fluctuations around the ordered classical 
phase can be quite significant for this system. Furthermore, we know that the intermolecular 
hopping amplitudes in the crystal are very sensitive to the interchain separation, and are 
somewhat sensitive the orientational state of the fullerene molecules. It would therefore be 
useful to carry out this analysis perturbing from the paramagnetic "metallic" phase of the 
doped system (containing both charge and spin fluctuations), as a function of varying the 
interchain separations. The structure of the resulting phase diagram as well as the character 
of the low lying quasiparticle excitations of the condensed phase in the presence of these 
ordered spiral spin textures would be particularly interesting to study. Even within the 
classical model, a distribution of J2/J1 ratios, such as might be introduced by orientational 
disorder on the fullerene sites, can be expected to couple strongly to the phase degree of 
freedom of the spiral phases and to lead to disordered phases with nontrivial magnetic 
properties which have yet to be studied. 

This work was supported in part by the Laboratory for Research on the Structure of 
Matter (University of Pennsylvania), by the NSF under the MRL program (Grant 92 20668) 

8 



and by the DOE (Grant 91ER 45118). Computations were carried out in part at the Cornell 
Theory Center, which receives major funding from NSF and New York State. 



9 



REFERENCES 

[1] O. Chauvet et ai, Phys. Rev. Lett. 72, 2721 (1994). 

[2] S.C. Erwin, G.V. Krishna and E.J. Mele, Phys. Rev. Lett, (submitted) (LU5465). 
[3] PW. Stephens et ai, Nature 370, 636 (1994). 

[4] W. Selke in Phase Transitions and Critical Phenomena, Vol. 15 (C. Domb and J. 
Lebowitz, ed; Academic Press, 1992). 



10 



FIGURES 

FIG. 1. (a) A portion of the pseudo-body-centered orthorhombic lattice for the polymerized 
phase, (b) The bco Brillouin zone and high-symmetry points in the q x — q z plane. 

FIG. 2. Spin structures on a cluster from the bco lattice. The axes of the polymer chains are 
indicated by the light lines. 

FIG. 3. Phase diagrams for the frustrated spin model on the bco lattice: (a) phases considering 
only colinear structures; (b) phase diagram for the Ising model; (c) phase diagram for the classical 
Heisenberg model. 

FIG. 4. Evolution of the pitch of the coplanar spiral phase as a function of the interaction ratio. 
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Fig. 1 E.J. Mele, G.V. Krishna, and S.C. Erwin 




Fig. 2 E.J. Mele, G.V. Krishna, and S.C. Erwin 




Fig. 3 E.J. Mele, G.V. Krishna, and S.C. Erwin 




Fig. 4 E.J. Mele, G.V. Krishna, and S.C. Erwin 



